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EXECUTIVE  SUMMARY 


This  is  one  of  five  companion  reports  prepared  under  the  sponsorship  of  Code  HECB  of  the  Air 
Force  Research  Laboratory  (originally  funded  by  Code  AL/OEBN,  Armstrong  Laboratory  of 
Wright-Patterson  Air  Force  Base).  Each  of  the  reports  deals  with  one  aspect  of  the  problem  of 
assessing  the  effects  of  noise  from  military  aircraft  on  marine  life:  (I)  criteria  and  thresholds  for 
injury  and  harassment.of  protected  marine  life,  (II)/(]II)  risks  of  impact  from  subsonic/supersonic 
aircraft  noise,  (IV)  metrics  for  sound  properties  in  air  compared  to  sound  properties  in  water,  and 
(V)  animal  population  statistics. 

The  end  purpose  of  this  multi-year  contract  effort  is  to  establish  technically  sound  estimation 
procedures  for  determining  the  effects  of  military  aircraft  noise  on  marine  life.  Without  such 
procedures,  the  Air  Force  risks  inadvertent  violations  of  the  law  and  becomes  vulnerable  to 
litigation  and  interference  with  military  operations. 

Objectives  of  the  contract  effort  include  developing  procedures  for: 

1)  Predicting  properties  of  sound  waves  in  air  and  under  water  as  generated  by  both 
subsonic  and  supersonic  aircraft  flights 

2)  Estimating  the  effects  of  sound  on  marine  life,  both  in  air  and  under  water 

3)  Determining  populations  of  marine  life  at  risk,  as  functions  of  aircraft,  flight  path,  and 
time  of  year. 

This  volume  specifically  addresses  the  approach  to  bounding  the  intensity  of  the  noise  field  in 
water  generated  by  a  supersonic  aircraft.  Note  here  that  the  propagation  mechanisms  for  aircraft 
noise  into  water  are  completely  different  for  sonic  booms  because  of  the  limited  source  angle; 
this  is  the  reason  for  a  separate  volume  on  the  subject. 

As  explained  in  Section  2,  the  properties  of  noise  in  water  of  interest  for  risk  assessment  are  the 
peak  pressure,  modified  positive  impulse,  and  energy  flux  density  in  selected  frequency  bands, 
for  each  boom  event.  The  number  of  events  over  time  is  also  important.  Multiple  exposure  rules 
vary  among  precedents  and,  just  as  for  impact  of  noise  on  humans  in  the  work  place,  may  not  be 
simple  energy  relationships. 

Section  3  reviews  the  well-studied  history  of  propagation  of  sonic  boom  noise  from  air  to  water, 
and  the  behavior  of  the  noise  in  water.  It  concludes  with  a  recommended,  approximate  approach 
for  readily  estimating  bounds  on  the  metrics  of  interest.  Practicality  is  emphasized,  especially  in 
dealing  with  the  rough  surface  interface  and  with  propagation  of  the  field  in  water. 

Section  4  provides  examples  for  different  aircraft  and  conditions.  It  shows  the  levels  of  the  noise 
in  deep  and  shallow  water.  Attention  is  also  paid  to  the  frequency  content  of  the  field,  since,  for 
example,  the  majority  of  marine  mammals  have  very  poor  hearing  at  low  frequency. 

The  most  important  conclusion  of  this  study  is  that  for  supersonic  flight  of  Air  Force  aircraft, 
there  are  very  few  and  limited  cases  for  which  there  might  be  a  compliance  issue  for  noise  in 
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water,  under  commonly  used  thresholds.  In  fact,  it  is  difficult  to  construct  examples  for  which 
there  is  a  statistical  risk  to  protected  marine  life,  even  for  focused  booms  and  exceptionally  good 
propagation  across  the  rough  air-sea  interface.  Such  cases  can  be  flagged,  and  care  taken  to 
avoid  areas  where  protected  species  may  be  present  and  at  risk  [A  companion  volume,  Moore 
and  Clarke  (1998),  provides  this  information  on  marine  species  at  a  summary  level  for  training 
areas,  including  most  of  the  US  coastal  region  to  200  nmi). 

The  principal  reason  for  the  lack  of  impact  of  sonic  boom  energy  under  water  is  that  even  for  the 
strongest  booms  and  good  coupling  to  the  water,  the  peak  pressure  and  energy  flux  density  are 
not  sufficient  to  cause  injury  or  harassment,  at  least  under  currently  accepted  criteria  and 
thresholds.  As  is  discussed  in  the  report,  an  extraordinarily  loud  boom  with  peak  pressure  as 
great  as  50  psf  at  the  water  surface  will  not  produce  the  12  psi  (1700  psf)  peak  pressure  in  water 
associated  with  harassment  of  marine  mammals  and  sea  turtles.  A  similar  situation  applies  for 
energy  thresholds. 

Nonetheless,  estimates  of  the  type  made  in  this  volume  are  essential  justifications  for  arguments 
of  insignificant  impact  (as  would  be  used  in  compliance  documents  or  regulator  consultations). 
They  also  anticipate  the  event  that  (as  happens  every  few  years)  the  favored  harassment  and 
injury  thresholds  are  changed  (from  the  view  of  the  regulators). 
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1.0  INTRODUCTION 


1.1  First  of  Series  of  Five  Reports 

This  is  the  third  in  a  series  of  five  companion  reports  prepared  under  the  sponsorship  of  the  Air 
Force  Research  Laboratory  (AFRL/HECB),  Wright-Patterson  Air  Force  Base  (originally  sponsored 
by  the  Noise  Effects  Branch,  Armstrong  Laboratory).  Each  of  the  reports  deals  with  one  aspect  of 
the  problem  of  assessing  the  effects  of  aircraft-generated  noise  on  marine  life: 

Report  I:  Criteria  and  Thresholds  for  Adverse  Effects  of  Underwater  Noise  on  Marine 
Animals 

Report  II:  Subsonic  Aircraft  Noise  at  and  beneath  the  Ocean  Surface:  Estimation  Models  for 
Metrics  Associated  with  Effects  on  Marine  Mammals 

Report  III:  Supersonic  Aircraft  Noise  at  and  Beneath  the  Ocean  Surface:  Estimation 
Models  for  Metrics  Associated  With  Effects  on  Marine  Mammals 

Report  IV:  Background  Definitions  and  Metrics  for  Sound  Properties  in  Air  and  in  Water 
Relevant  to  Noise  Effects 

Report  V:  Marine  Animal  Populations  for  Ocean  Regions  of  Interest  to  Air  Force  Flight 
Operations 

The  end  purpose  of  this  multi-year  contract  effort  is  to  establish  technically  sound  estimation 
procedures  for  determining  the  effects  of  military  aircraft  noise  on  marine  life.  Without  such 
procedures,  the  Air  Force  risks  inadvertent  violations  of  the  law  and  becomes  vulnerable  to 
litigation  and  interference  with  military  operations. 

Objectives  of  the  contract  effort  include  developing  procedures  for: 

1)  Predicting  properties  of  sound  waves  in  air  and  under  water  as  generated  by  both 
subsonic  and  supersonic  aircraft  flights 

2)  Estimating  the  effects  of  sound  on  marine  life,  especially  under  water 

3)  Determining  populations  of  marine  life  at  risk,  as  functions  of  aircraft,  flight  path,  and 
time  of  year. 

This  volume  deals  with  supersonic  aircraft  noise. 


1.2  Organization  of  Report  and  Appendices 

Section  2  of  this  report  is  intended  to  provide  context  and  motivation  for  the  effort.  In  particular, 
it  gives  examples  of  the  types  of  information  about  the  noise  field  that  are  needed  to  assess  risk 
to  marine  life.  Perhaps  the  most  important  property  common  to  injury  and  harassment  thresholds 
for  protected  species  is  the  dependence  on  exposure  time  (number  of  intermittent  events). 
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Section  3  contains  a  review  of  the  literature  on  the  subject  of  sound  propagation  from  a  source  in 
air  to  a  receiver  at  or  below  the  air-sea  interface.  From  this  review,  practical  approaches  are  then 
recommended  for  estimating  the  acoustic  field  from  an  aircraft  traveling  at  supersonic  speeds. 
Note  that  noise  from  subsonic  sources  has  different  propagation  physics  and  is  covered  in  a 
separate  report. 

Section  4  gives  examples  of  the  time  series  of  sound  pressure  level  that  would  be  observed  at  a 
point  in  the  ocean,  as  generated  by  a  variety  of  aircraft  types  and  conditions. 

An  important  conclusion  is  drawn:  it  would  be  a  rare  event  for  a  supersonic  aircraft  to  generate 
sufficient  noise  in  water  to  be  considered  a  source  of  harassment  of  protected  species,  at  least  for 
commonly  used  criteria  and  thresholds.  Furthermore,  the  approach  of  this  volume,  combined 
with  data  from  companion  volumes  on  criteria/thresholds  and  on  mammal/turtle  populations,  can 
be  used  to  identify  those  cases  of  concern  and  afford  Air  Force  the  opportunity  to  greatly  limit 
risk. 


1.3  Acknowledgments 
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2.0  GEOMETRIES  AND  METRICS  OF  INTEREST 


The  purpose  of  this  report  is  not  to  advance  modeling  or  data  analysis  technology  of  aircraft 
noise  in  water,  but  rather  quite  specifically  to  provide  approaches  and  examples  of  how  to  predict 
the  properties  of  aircraft  noise  in  water  needed  to  assess  compliance  risk.  The  types  of  noise 
properties  needed  for  this  problem  (risk  assessment  for  impulsive  noise)  include  at  least  three 
different  quantities:  peak  pressure,  energy  flux  density,  and  modified  positive  impulse,  with  rms 
pressure  an  additional  possibility.  Even  these  metrics  can  have  complicated  constraints,  such  as 
“the  largest  1/3 -octave-band  energy  level  above  100  Hz”  from  a  recent  Navy  EIS.  Some 
examples  are  provided  below,  to  help  identify  what  information  may  be  needed  and  the  scales  of 
the  problem. 

Additionally,  geometric  relationships  among  the  aircraft,  the  sea  surface,  and  the  marine  animals 
are  outlined  for  context. 


2.1  Properties  of  the  Noise  Field  Needed  for  Risk  Assessment 

A  companion  volume  to  this  report  summarizes  criteria  and  thresholds  for  injury  and  harassment 
of  protected  marine  species  by  underwater  noise.  Precedents  in  the  eyes  of  the  regulators  and 
other  DOD  branches  (especially  Navy  and  DARPA)  must  be  acknowledged,  along  with  the  view 
of  the  scientific  community. 

For  impulsive,  underwater  noise  (as  for  explosives,  airguns,  and  sonic  booms),  injury  and 
mortality  thresholds  for  marine  mammals  and  sea  turtles  currently  used  (and  likely  to  be 
applicable  in  the  future)  are  of  form: 

•  peak  pressure  (including  a  penalty  for  multiple  exposures), 

•  energy  flux  density  in  a  frequency  band  (also  depending  on  exposure  count), 

•  modified  positive  impulse  (with  modification  depending  on  animal  depth  in  water) 

Bandwidth  modifiers  generally  reflect  the  hearing  band  of  the  species  in  question,  just  as  A- 
weighting  is  used  in  evaluating  noise  impact  on  humans.  There  is  also  the  complication  of  special 
effects  thought  to  be  related  to  low-frequency  sound. 

The  table  below  provides  some  relevant  examples  of  thresholds  used  in  recent  compliance  work. 
More  on  this  topic  can  be  found  in  the  companion  volume  (Cavanagh  and  Laney,  2000).  In  the 
table,  decibel  quantities  for  energy  flux  density  are  referenced  to  1  |iPa2-s. 
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Table  2-1.  Examples  of  Thresholds 


Harassment  Criterion 

Threshold 

Reference 

TTS  for  all  small  odontocetes 

Energy  flux  density  level  for 
band  above  100  Hz  >  182  dB  - 
17  log  N,  where  N  is  number  of 
impulses.  Lower  threshold  limit 
=120  dB. 

SEAWOLF  Shock  Trial 

FEIS  (1998)  and  NMFS 

Final  Rule  (1998) 

Injury  for  sea  turtles 

Peak  pressure  above  15  psi 

Young  (1991) 

Injury  for  mammals  in  the 
form  of  eardrum  rupture 

Energy  flux  density  in  excess  of 
1.2  psi-in  (205  dB) 

SEAWOLF  Shock  Trial 

FEIS  (1998) 

PTS  for  marine  mammals 

rms  pressure  level  exceeds  190 
dB  (re  1  (iPa) 

HESS  committee,  as 
discussed  at  NMFS  criteria 
workshop  (1998) 

Here  TTS  and  PTS  are  degradations  in  hearing,  sometimes  treated  by  regulators  as  harassment 
criteria  for  mammals  under  the  Marine  Mammal  Protection  Act  and  Endangered  Species  Act. 

To  avoid  risk  or  to  show  that  there  is  no  impact,  estimates  are  needed  of  the  relevant  properties 
of  noise  as  received  by  the  animal.  Peak  pressure  and  energy  in  a  band  are  the  most  common 
metrics.  This  paper  then  provides  approaches  for  estimating  these  exposure  values  in  the  water 
for  the  case  of  supersonic  aircraft  sources. 


2.2  Animal/Interface/Aircraft  Geometries 

In  assessing  risk  of  injury  or  harassment  of  protected  species,  the  key  elements  are  the  time 
history  of  the  noise  at  a  site  in  the  ocean  (as  outlined  above)  and  the  movement  of  the  animals  of 
concern.  The  former  is  a  function  of  aircraft  source  level,  speed,  altitude,  and  lateral  distance,  as 
well  as  of  the  propagation  environment  in  air  and  in  water.  At  each  point  in  the  ocean,  the  time 
series  can  be  estimated,  and  the  impact  assessed.  Statistical  descriptions  of  animal  populations 
are  then  applied  to  determine  the  likelihood  of  encounter. 

For  multiple  exposure  potentials,  as  perhaps  from  formation  flights  or  multiple  explosives, 
animal  motion  most  also  be  addressed.  The  problem  often  comes  down  to  one  of  determining 
the  chances  that  an  animal  will  stay  within  an  area  for  a  specific  amount  of  time.  Rather  than 
address  the  animal  motion  problem  here  (it  takes  many  pages),  we  can  only  alert  the  planner  to 
the  problem,  and  note  that  several  methods  of  solution  have  been  used  in  the  past  (see,  for 
example,  SEAWOLF  (1998)  and  SSQ-1 10  EA  (1995)). 

The  usual  approach  to  risk  estimation  for  moving  sources  and  exposure-time-dependent 
thresholds  for  injury/harassment  is  to  map  out  the  ocean  region  consisting  of  all  locations 
ensonified  at  or  above  threshold  levels.  The  expected  number  of  animals  affected  is  then  the  size 
of  that  ensonified  region  multiplied  by  animal  density.  The  ensonified  region  is  usually 
characterized  in  terms  of  area,  and  terms  such  as  ‘footprint’  or  ‘sweep  region’  are  sometimes 
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used.  If  there  is  a  depth  bias  for  the  ensonification  and  the  animal  location,  then  the  appropriate 
volume  of  the  region  and  revised  animal  densities  must  be  considered. 

As  an  illustrative  example,  suppose  an  aircraft  ensonifies  the  ocean  region  within  1000  m  on 
either  side,  to  depths  of  1000  m,  at  threshold  levels  (e.g.,  energy  level  in  excess  of  182  dB).  If 
the  aircraft  follows  a  non-overlapping  path  for  one  hour  at  200  km/hr,  then  the  impact  region 
would  have  area  of  200  km2.  If  the  density  of  blue  whales  was  0.0004/km2 ,  then  the  expected 
number  of  blue  whales  to  be  affected  would  be  about  0.08,  well  below  what  would  ordinarily  be 
considered  a  significant  risk. 


2-3 


5 


3.0  NOISE  ESTIMATION 


This  section  reviews  the  standard  references  for  sonic  boom  generation  and  propagation, 
and  then  estimates  corresponding  noise  properties  at  the  air-sea  interface  and  in  the  water 
column.  Noise  properties  of  special  interest  are  peak  pressure,  pulse  duration, 

“footprint,”  energy  spectrum,  and  dependencies  on  aircraft  kinematics  and  receiver  depth. 

3.1  Review  of  Features  of  Sonic  Boom  Generation  and  Propagation  in  Air 

Objects  moving  in  the  atmosphere  at  speeds  greater  than  the  local  speed  of  sound 
generate  shock  waves  or  sonic  booms.  In  the  case  of  a  simple  projectile,  two  primary 
waves  are  generated,  one  at  the  front  of  the  object  and  the  other  at  the  rear  (Figure  3-1). 
Aircraft  present  a  more  complicated  picture  with  waves  coming  from  the  bow,  air 
intakes,  control  surfaces,  and  tail.  At  large  distances  from  the  aircraft,  however,  the 
system  of  shock  waves  tends  to  distort  and  steepen  until  they  coalesce  into  a  bow  and  tail 
wave  similar  to  that  from  the  simple  projectile  (Figure  3-2). 


Figure  3-1  Far-field  Sonic  Boom  Pattern  (Maglieri  &  Plotkin,  1995) 


Figure  3-2  Boom  Signature  as  Function  of  Altitude  (Maglieri  &  Plotkin,  1995) 
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As  the  bow  wave  approaches,  the  ambient  pressure  rises  quickly  to  Ap  above  atmospheric 
pressure.  Following  the  bow  wave  a  decompression  occurs  which  is  corrected  with  a 
sudden  recompression  from  the  tail  wave.  Generally  the  bow  and  tail  waves  have  similar 
strength  and  pressure  decreases  linearly  between  the  two  waves.  This  sudden 
compression,  linear  decompression,  and  sudden  recompression  form  the  classic  sonic 
boom  N-wave.  N-waves  sweep  out  a  sonic  boom  carpet  whose  width  depends  on  flight 
and  atmospheric  conditions  (Figure  3-3).  The  primary  boom  is  confined  to  a  conical 
region  extending  behind  the  aircraft.  The  primary  boom  footprint  on  the  ground  is  a 
hyperbolic  intersection  of  the  boom  cone  with  the  ground  (Figure  3-4).  The  primary 
boom  cone  forms  an  angle  of  incidence,  a,  measured  from  the  normal  vector  to  the  air- 
water  interface  computed  as  follows: 


or  =  sm 


-i 


XMJ 


where  M  is  the  aircraft  Mach  number  and  M>1  (Figure  3-5).  As  described  in  Section  3.3, 
this  angle,  known  as  the  Mach  angle  has  value  90°  for  M=l,  30°  for  M=2  and  approaches 
0°  as  M— For  level  flight  and  M<4.35,  a  is  greater  than  the  critical  angle  (13.3°)  for 
passage  of  sound  from  the  air  to  the  water  and  is  the  most  important  factor  in  determining 
how  sonic  boom  energy  couples  into  the  ocean. 


Figure  3-3  Nature  of  Sonic  Boom  Ground  Exposure  (Maglieri  &  Plotkin,  1995) 

Following  a  supersonic  flyby,  the  observer  on  the  ground  may  hear  two  booms  in  quick 
succession.  The  first  boom  corresponds  to  the  bow  wave  and  the  second  to  the  tail  wave. 
For  small  objects,  the  two  booms  may  merge  into  a  single  detectable  event.  For  aircraft, 
the  time  interval  between  the  bow  and  tail  wave  is  dependent  on  the  length  of  the 
airframe,  aircraft  altitude,  and  aircraft  speed.  For  most  aircraft  the  time  interval  is  0.1 
seconds  or  more,  allowing  the  observer  to  detect  two  distinct  booms. 
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Figure  3-4  Sonic  Boom  Cone  (Carlson  &  Maglieri,  1972) 


Figure  3-5  Wave  and  Mach  Cone  Geometry 


In  addition  to  the  primary  boom  carpet  characterized  by  the  direct  path  N-wave,  a 
secondary  boom  carpet  may  exist.  The  extent  of  the  primary  boom  carpet  is  limited  by 
the  lateral  cutoff  imposed  by  the  sound  speed  profile  of  the  atmosphere.  Figure  3-6 
shows  typical  wind  and  temperature  profiles  for  the  atmosphere.  These  profiles  limit  the 
range  of  direct  downward  path  propagation  from  a  source  at  some  altitude  to  the  ground. 
The  secondary  boom  carpet  is  formed  by  waves  which  either  bounce  in  the  region  of  the 
primary  carpet  and  are  then  bent  down  in  the  upper  atmosphere  or  by  waves  which 
propagate  upward  and  are  bent  down  in  the  upper  atmosphere  into  the  region  of  the 
secondary  carpet  (Figure  3-6).  Secondary  booms  are  weak  in  peak  pressure  (0.02  -  0.20 
lb/ft2  compared  to  0.5-  3.5  lb/ft2  for  typical  primary  booms)  and  are  in  the  0.1  to  2.0  Hz 
range. 


Variations  in  the  atmosphere  also  affect  the  sonic  boom  signature  shape.  Turbulence  and 
thermal  variations  in  the  lower  atmosphere  often  distort  the  sonic  boom  waveform 
leading  to  modified  N-waves  which  are  peaked  or  rounded  at  either  the  leading  or  trailing 
edge  of  the  N-wave  or  at  both  ends  (Figure  3-7).  The  categories  in  Figure  3-7  were 
developed  by  NASA  and  the  labels  refer  to  the  extent  and  position  of  observed 
“peaking”,  “spiking”,  and  “rounding”.  When  peaking  or  spiking  occurs,  maximum 
overpressures  higher  than  predicted  for  a  regular  N-wave  are  observed. 


Figure  3-7  Categories  of  Sonic  Boom  Waveforms  (Pierce,  1968) 


3.2  Noise  Properties  at  the  Air-Sea  Interface 

Numerous  studies  have  been  conducted  to  predict  and  measure  the  pressure  signatures 
associated  with  the  primary  boom  carpets  produced  by  objects  in  supersonic,  steady,  level 
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flight.  For  aircraft,  the  measurements  and  theory  agree  with  boom  levels  increasing  with 
aircraft  size  and  decreasing  with  increasing  altitude  (Figure  3-8).  Peak  overpressures  for 
sonic  booms  are  generally  in  the  0.5-  3.5  lb/ft2  range  and  are  traditionally  expressed  in 
lb/ft  in  the  literature.  For  reference,  a  peak  pressure  of  1  lb/ft2  is  approximately 
4.758-107  (iPa  =  47.58  Pa  =  6.94-10'3  psi  =  153.5  dB  (re  1  pPa)  =  127.5  dB  (re  20pPa). 


6r 


O  Aircraft  data 
#  Ascent  launch  vehicles 
■  Apollo  capsule  reentry 
A  Shuttle  orbiter  reentry 
Prediction 


Sonic  boom  4  [~  XB-70 
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Figure  3-8  Measured  and  Predicted  Overpressures  (Maglieri  &  Plotkin,  1995) 


Assuming  steady  level  flight  and  N-wave  arrivals,  the  following  approximations  from 
Carlson,  1978  are  generally  correct  to  within  5%  of  computer  calculations  (Maglieri  & 
Plotkin,  1995).  Carlson  also  includes  approximations  for  off  ground  track,  climbing,  and 
descending  peak  pressures  which  are  not  presented  here.  (Note  that  in  Carlson’s  paper, 
the  units  for  effective  aircraft  altitude  he  should  be  in  meters  for  computing  Apmax  and  At.) 


=1K„K,JWs{M2-\)'‘hf'l 

3.42  M 


3/4 


A t  =  K. 


Pv~Pge 


c 
-0.12 


(M2-lf 


h!,4i3'*K . 


/8  "e 


Wm2-  i 

1.4PVM2/2 


[1] 


where 

Ap 

max  shock  strength  (Pa) 

At  N-wave  duration  (s) 

Pv  ambient  pressure  at  vehicle  altitude  (Pa) 

Pg  ambient  pressure  at  ground  (Pa,  nominally  100  kPa ) 

cv  speed  of  sound  at  aircraft  altitude  (nominally  300  m/s) 

h  aircraft  altitude  above  sea  level  (km) 

he  effective  aircraft  altitude  above  local  ground  level  (m) 
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Figure  3-10  Signal  Duration  Factor  (Carlson,  1978) 


The  power  spectrum  for  an  N-wave  is  computed  as  follows  (Howes,  1967): 


\p{p\ 2  =  ~^4£)2~j  I4  “  4  cos  (cqD) + ( coD f  [l  +  cos(tyZ))]  -  4  coD  sin(<y£>)} 


where  to  is  sinusoidal  frequency  and  D  is  the  N-wave  duration  in  seconds. 

Asymptotically  the  power  spectral  density  changes  by  6  dB/octave  with  positive  slope  in 
the  low  frequencies  and  negative  slope  in  the  high  frequencies  around  the  peak  (Figure  3- 
11).  It  should  be  noted  that  the  energy  spectral  density  equation  presented  by  Young, 
1966  is  incorrect.  The  frequency  of  the  first  and  dominant  peak  in  the  power  spectrum  is 
given  by  Howes,  1967  as: 

fmax  ~  2/(3D). 


Since  typical  N-wave  durations  are  0.035  -  0.3  seconds,  typical  fmax  are  in  the  2  -  20  Hz 
range.  Assuming  an  8  Hz  peak,  the  spectrum  level  at  1000  Hz  will  be  over  40  dB  down- 
from  the  peak  spectrum  level. 


Figure  3-11  First  Peak  in  Power  Spectra  for  Various  Aircraft  (Young,  1966) 

In  addition  to  the  peak  magnitude  and  frequency  content  of  a  sonic  boom,  its  lateral 
spread  is  also  of  interest.  “The  extent  of  the  primary  carpet  is  the  point  at  which  the  ray 
refracts  away  from  the  ground  (cutoff  distance)”  (Magielri  &  Plotkin,  1995).  This  cutoff 
point  is  determined  by  the  characteristics  of  the  atmosphere,  aircraft  altitude,  and  aircraft 
speed.  Figure  3-12  shows  boom  overpressures  as  a  function  of  lateral  distance  from  the 
ground  track  for  an  aircraft  at  an  altitude  of  60,000  ft  traveling  in  steady  level  flight  at 
approximately  M=2.0.  “The  widths  of  the  sonic  boom  carpets  on  the  ground  increase 
with  increasing  altitude  and  Mach  number”  (Magielri  &  Plotkin,  1995).  Figure  3-13 
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compares  the  theoretical  and  measured  lateral  extent  of  sonic  booms  as  a  function  of 
aircraft  speed  and  altitude.  “A  crude  but  useful  rule  of  thumb  equates  the  lateral  spread  in 
miles  to  the  airplanes  altitude  in  thousands  of  feet”  (Carlson  &  Maglieri,  1972). 


Signature  shape 


Figure  3-12  Overpressures  as  a  Function  of  Lateral  Distance,  M=2  (Maglieri  & 

Plotkin,  1995) 
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Figure  3-13  Sonic  Boom  Width  on  Ground  (Maglieri  &  Plotkin,  1995) 


It  is  useful  to  summarize  what  the  observer  on  the  ground  would  experience  as  the  result 
of  a  steady  level  supersonic  flyby.  An  observer  directly  below  the  flight  path  would 
experience  an  N-wave,  possibly  distorted  as  in  Figure  3-7,  with  a  peak  pressure  and 
duration  which  can  be  estimated  using  equations  [1].  As  an  observer  moves  laterally 
away  from  the  ground  track,  the  signature  duration  will  increase  and  the  peak  pressures 
will  decrease  as  in  Figure  3-12.  As  the  observer  moves  further  from  the  ground  track,  the 
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signal  will  disappear  as  the  observer  crosses  the  lateral  cut-off  defining  the  primary  boom 
carpet  (Figure  3-3).  Moving  further,  an  observer  would  experience  the  secondary  boom 
carpet  where  peak  pressures  are  less  than  10%  of  the  primary  boom  peak  pressures  and 
where  energy  is  below  several  Hertz.  Typical  parameters  for  a  ground  observer 
following  a  steady  level  flight  are  shown  in  Table  3-1 . 


Parameter 

Primary  Boom 

Secondary  Boom 

Peak  Pressure  (lb/ft2) 

0.5 -3.5 

0.02-0.2 

Peak  Pressure  (dB  re  20pPa) 

121.5-138.4 

93.5-113.5 

Wave  Duration  (sec) 

0.035  -  0.3 

>  10.0 

Wave  Shape 

N-wave 

variable 

Boom  Carpet  Width  (nmi) 

20-60 

20  -  50 

Mach  Angle,  a  (degrees) 

20-65 

N/A 

Peak  Frequency  (Hz) 

2-20.0 

0.1  -2.0 

Table  3-1  Typical  Boom  Parameter  Values 


As  with  processes  such  as  weathering  and  corrosion,  sonic  booms  do  not  generally  create 
instantaneous  damage  to  buildings  and  other  structures  but  rather  contribute  to 
deterioration  over  time  (Warren,  1972).  Table  3-2  describes  the  effect  on  structures  and 
public  reaction  to  sonic  booms  with  various  levels  of  peak  overpressure  (Newman  and 
Beattie,  1985). 


Peak  Overpressure  (psf) 

Damage  to  Ground  Structures 

Public  Reaction 

0-1 

None 

None 

|  1-11-5 

None 

Probable 

1.5-1.75 

None 

Significant 
(especially  at  night) 

1.75-2.0 

None 

Significant 

2.0-3.0 

Incipient 

Widespread 

Table  3-2  Typical  Sonic  Boom  Damage  and  Public  Reaction  (Newman  and  Beattie, 

1985) 


3.3  Noise  Properties  Across  the  Air-Sea  Interface  and  in  the  Water 

This  section  addresses  the  problem  of  propagation  of  a  sonic  boom  wave  from  the  air 
across  the  air-sea  interface  into  the  ocean.  Treatments  in  the  literature  sometimes  appear 
to  be  overly  complicated,  with  terms  such  as  the  ratio  of  aircraft  speed  to  sound  speed  in 
water  (Mach  number  in  water). 

In  this  review,  the  problem  is  considered  in  two  parts: 
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(a)  propagation  of  harmonic  plane  waves  from  air  to  water,  the  critical  angle,  and 
the  special  case  of  plane  waves  arriving  at  angles  shallower  than  critical  angle, 
and 

(b)  properties  of  the  sonic  boom  wave  in  air  at  the  interface  (especially  its 
amplitude,  waveform,  and  arrival  angle  (Mach  angle)),  and  application  of  (a)  to 
find  the  refracted  field  in  water. 

Topic  (a)  is  straightforward,  and  the  approach  applies  to  harmonic  plane  waves  in 
general.  The  results  have  been  known  since  at  least  the  time  of  Rayleigh  (Strutt,  1848), 
and  are  well  explained  in  the  classic  textbooks  of  Officer  (1958)  and  Brekovskikh  (1960). 
The  physical  principles  are  quite  general,  and  do  not  distinguish  between  sonic  booms 
and  other  plane  waves. 

Topic  (b),  on  the  other  hand,  uses  the  results  of  (a)  and  fourier  synthesis  to  estimate  the 
N-wave  field  in  water.  The  fourier  method  is  routinely  used  in  underwater  sound  to 
estimate  the  propagation  of  a  finite-bandwidth  waveform  from  the  properties  of  the 
propagation  of  pure  tones,  which  properties  are  readily  available  from  popular  wave 
models  such  as  normal  mode  and  parabolic  equation  models.  The  fourier  approach  is 
older  than  Rayleigh,  and  useful  explanations  for  the  sonic  boom  case  are  given  in 
Sawyers  (1968)  and  Cook  (1970).  Note  that  the  waveform,  amplitude,  and  arrival  angle 
(Mach  angle)  are  the  only  properties  of  the  sonic  boom  used.  The  Mach  number  is 
relevant  only  in  its  determination  of  those  properties,  and  need  not  appear  at  all  in  the 
formulas  for  the  field  in  water. 

As  noted  often,  the  limited  impact  of  sonic  boom  noise  in  water  is  a  result  of  the  special 
coincidence  that  the  Mach  angles  for  most  aircraft  speeds  are  not  steep  enough  for  good 
transmission  into  water.  If  water  had  a  lower  sound  speed  or  airplanes  traveled  at  Mach 
5,  there  would  be  a  significantly  greater  impact  of  the  sonic  boom  in  water. 


3.3.1  Harmonic  Plane  Wave  Reflection  and  Transmission 

From  the  beginning  pages  of  Brekhovskikh  (1960,  pp  16  ff),  consider  the  Rayleigh 
reflection  relationships: 

Pinc  =  A  exp  (ika(xsin0a  -  zcos0a)) 
pref  =  VA  (exp  (ika(xsin0a  +  zcos0a)) 

Pw  =  Ptr  =  WA  exp  (ikw(xsin0w  -  zcos0w)) 


where 


p  refer  to  the  incident,  reflected  and  transmitted  pressures, 
subscript  a  is  for  air  and  w  is  for  water. 


3-10 

15 


c  is  sound  speed,  0  is  angle  (measured  from  the  vertical), 
(O  is  angular  frequency 
k  =  co/c  is  wavenumber 

A  =  -ipco  (amp),  where  amp  is  the  pressure  amplitude 
and 

the  frequency  factor  exp(icot)  has  been  omitted. 


From  the  boundary  conditions  of  continuity  of  pressure  and  continuity  of  the  normal 
component  of  the  particle  velocity  across  the  interface,  the  Rayleigh  coefficients  are 
found  as: 


V  =  [m  cos  0a  -  (n2  -  sin20a)1/2  ]/[m  cos  0a  -  (n2  -  sin20a)1/2] 


W  =  (l/m)(l+V) 


where 


m  =  pw/pa  and  n  =  ca/cw  =  kw/ka 


3.3.2  Plane  Wave  Propagation  form  Air  to  Water,  and  the  Critical  Angle. 

For  propagation  from  air  into  water,  m  has  value  of  about  900  and  n  about  0.23.  Hence, 
for  incidence  angles  0a>  arcsin(n)  =  13  degrees,  V  is  a  complex  number  with  value  close 
to  one  (perfect  reflection).  However,  the  transmitted  pressure  in  water  satisfies: 

Pw  =  Ptr  =  WA  exp  [ikw(xsin0w  -  zcos0w)] 

=  (l/m)(l+V)A  exp[ikaxsin0a  +  kaz(n2  -  sin20a)1/2] 

This  corresponds  to  a  wave  traveling  in  the  x  direction  with  exponential  decay  in  the  z 
direction.  In  fact, 

I  pw  I  =  (2/m)  |  A  |  exp[-kaz  (n2  -  sin20a)1/2] 

The  angle  arcsin(n)  =13  degrees  is  the  critical  angle,  so  that  any  plane  wave  in  air  that 
arrives  at  the  air-sea  interface  at  an  angle  less  steep  than  about  13  degrees  will  not  be 
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transmitted  into  the  water  as  a  propagating  wave,  but  rather  as  an  evanescent  wave,  i.e., 
one  which  decays  exponentially  with  depth  in  the  water. 


WATER  PROPAGATING  EVANESCENT 

ARRIVALS  arrivals 

Figure  3-14  Critical  Angle  Geometry 


3.3.3  Application  to  Sonic  Booms  -  Mach  Angle  vs  Critical  Angle 

As  discussed  in  section  3.1),  an  aircraft  in  level  flight  at  speed  Mach  M  >  1  (i.e.,  at  speed 
Mca )  produces  a  shock  wave  at  the  Mach  angle  (a  =  arcsin(l/M)).  From  Figure  3-5,  note 
that  this  is  the  angle  measured  from  the  axis  of  the  aircraft  path,  or  equivalently  the  angle 
that  the  shock  wave  makes  with  the  vertical.  It  is  thus  the  incidence  angle  of  the  shock 
wave  at  the  air-sea  interface,  and  can  be  compared  directly  with  the  critical  angle  of 
3.3.2. 

Some  sample  values  of  a  illustrate  that  for  1<M<4.3,  the  incidence  angle  is  less  steep 
than  the  critical  angle,  and  the  wave  filed  will  be  transmitted  into  the  water  only  as  an 
evanescent  wave: 


3.3.4  Application  to  Sonic  Booms  -  Synthesized  N-Waves 

To  finally  address  sonic  boom  propagation,  begin  with  an  N-wave  incident  on  the  air-sea 
surface  at  angle  a.  The  results  of  subsection  3.3.2  apply  to  harmonic  waves  (pure  tones), 
but  give  the  pressure  field  below  the  interface  as  a  function  of  arrival  angle,  position  in 
the  water,  and  frequency. 

For  aircraft  speeds  less  than  about  M=4.3,  the  Mach  angle  a  will  be  greater  than  the 
critical  angle.  Thus,  the  incident  sonic  boom  wave  in  the  air  is  reflected  from  the  surface. 
The  only  energy  entering  the  water  is  in  the  form  of  subsonic,  non-radiating  plane  waves 
whose  amplitudes  ‘evanesce’  or  decay  exponentially  with  depth  below  the  surface.  Air 
Force  manned  aircraft  fly  below  M=4  and  so  all  associated  sonic  booms  have  only  an 
evanescent  transmission  into  the  water. 


For  an  aircraft  in  constant  supersonic  level  flight  over  a  flat  air-water  interface,  the 
underwater  pressure  field  generated  by  a  theoretical  perfect  N-wave  can  be  computed  via 
the  fourier-synthesis  method  mentioned  earlier.  Following  Cook,  1970:  “The  first  step  is 
to  find  the  reflected  and  refracted  waves  for  an  incident  sinusoidal  wave.  The  second 
step  is  to  find  the  Fourier  transform  of  the  incident  N-wave,  and  so  to  determine  the  total 
reflected  and  refracted  waves  (caused  by  the  N-wave)  by  linear  superposition  of  the 
effects  of  the  incident  sinusoidal  components”.  Following  this  prescription  leads  to  the 
following  underwater  sound  pressure  estimate  for  an  incident  N-wave  where  M  <  W: 


where 


pw(x',h)  =  2p0  cos  A(/j  cos  A + 1 2  sin  A) 
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s  =  {p/p.)(i-M2 /w2f2  /{m1  -if 


-2 


12 


[2] 


x  =  scaled  horizontal  distance  from  center  of  N-wave  (units  of  lh  N-wave  length) 
z  =  scaled  depth  below  surface  (units  of  Vi  N-wave  length) 

M  =  aircraft  Mach  number 

W  =  ratio  of  speed  of  sound  in  water  to  that  in  atmosphere  (—4.5) 
p  =  density  of  air  (-1.21  kg/m3) 
pw  =  density  of  water  (-1000  kg/m3) 
po  =  peak  pressure  of  incident  N-wave 
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The  (x,y,z)  components  of  the  wavenumber  vector  for  an  incident  sinusoidal  wave  are 
(ko,0,kz).  The  magnitude  of  this  vector  is  k=<n/c  where  c  is  the  speed  of  sound  in  air.  The 
same  components  for  the  underwater  inhomogeneous  refracted  wave  are  (ko,  0,  -iK).  The 
complete  set  of  wave  numbers  (in  terms  of  k=co/c)  for  an  incident  sinusoidal  wave  is  as 
follows: 

k0  -klM  =  ksma 
kz=(l-l/M2f2k  =  kcosa 
K  =  (l/M2-\/W2f2k 

In  terms  of  (x’,y’,z’)  the  wavenumbers  become: 

k0  =  k 

kz=(M2-lf2k 
K  =  (\-M2/W2f2k 

where  K/k  =  (I  from  equations  [2] .  The  amplitude  ratio  of  the  refracted  to  incident 
sinusoidal  wave  for  M<W  can  be  written  as  (Officer,  1958): 


where  A  from  equation  [2]  is  seen  to  be  the  phase  change  between  the  incident  and 
inhomogeneous  refracted  wave.  Also,  the  evanescent  decay  of  a  sonic  boom  underwater 
scales  as  e~kofK’  (as  shown  in  section  3.3.2). 

For  most  practical  situations  cosA~1.0  and  sinA=0  so  pw=2Iipo.  For  example,  if  M=1.5 
then  8=0.001,  A=0.06,  cosA=  1.000  and  sinA~  001.  Given  this  result,  the  main  part  of  the 
sound  field  is  given  by  Ii(x’,h).  This  function  is  an  N-wave  at  the  surface  and  spreads  out 
underwater  to  form  infinitely  long  precursor  and  tail  waves  (Figure  3-15).  I2(x’,h)  has 
infinite  spikes  at  x’=  ±1,  h=0  and  also  has  infinitely  long  tails  and  precursors.  I2(x’,h) 
however  makes  a  negligible  contribution  to  the  total  sound  pressure  beyond  a  few 
centimeters  from  x’=  ±1  for  most  cases  (Cook,  1970). 

Sparrow  (1995)  points  out  that  equations  [2]  show  that  “faster  aircraft  speeds  will 
produce  higher  peak  pressures  at  a  fixed  depth”  (all  other  things  being  equal,  the  greater 
the  speed,  the  steeper  the  sonic  boom  arrival).  Figure  3-15  shows  the  normalized 
pressure  resulting  from  an  N-wave  on  the  ocean  surface  as  a  function  of  distance  from  the 
center  of  the  N-wave  and  depth  as  predicted  by  equation  [2]  where  M=2.7  .  Note  that  at 
the  ocean  surface  the  peak  pressure  is  twice  that  of  the  incident  wave  in  air  and  that  the 
peak  pressure  drops  by  approximately  a  factor  of  8  at  a  depth  of  half  the  N-wave  length. 
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The  N-wave  travels  at  approximately  M*333m/s  where  M  is  the  aircraft  Mach  number. 
For  typical  aircraft,  N-wave  durations  are  0.1  -  0.3  seconds  leading  to  underwater 
exposure  times  of  less  that  one  second. 


Scaled  Presure  of  Incident  N-wave 


Figure  3-15  Example  of  Underwater  Sound  Pressure  as  Estimated  by  Equations  [2] 


As  an  N-wave  penetrates  the  ocean,  high  frequency  components  are  quickly  attenuated. 
An  example  is  provided  below  (Figure  3-16)  for  an  aircraft  traveling  Mach  2.2  creating  a 
2.5  psf  peak  pressure  N-wave  with  a  duration  of  0.3  seconds.  Note  that  by  a  depth  of  4.6 
m,  that  N-wave  pressure  levels  fall  below  typical  ambient  noise  levels  below  200  Hz.  It 
should  be  noted  that  the  “typical”  ambient  levels  presented  in  Figure  3-16  are  slightly 
higher  than  the  values  reported  by  Urick,  1983.  Also,  the  levels  presented  in  Figure  3-16 
are  spectrum  levels. 
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PENETRATION  OF  A  SONIC  BOOM  INTO  THE  OCEAN 
(adapted  from  Cook  it  al.,  1972) 


Figure  3-16  Sound  Spectra  of  Sonic  Boom  at  and  Below  Ocean  Surface  (Cummings, 

1994) 


3.4  Special  Cases:  Aircraft  Maneuvers,  Focused  Beams,  Rough  Surface 

In  the  case  of  a  diving  aircraft,  an  N-wave  arriving  at  less  than  the  critical  angle,0,  can  be 
produced  (see  Figure  3-14  for  geometry).  For  plane  waves  incident  at  angles  less  than 
critical,  the  pressure  in  the  water  is  approximately  twice  the  pressure  of  the  wave  in  the 
air.  The  transmitted  intensity,  however,  is  only  about  0.1 1%  of  the  incident  intensity  due 
to  the  impedance  mismatch  between  the  air  and  water.  This  corresponds  to  an  intensity 
decrease  of  approximately  30  dB.  A  plane  wave  with  an  incident  angle  0a,  less  than 
critical,  refracts  into  a  plane  wave  in  the  water  with  initial  incidence  angle  0w  related  by 
sin0w  =  (cw/ca)  sin0A  where  cw>a  is  the  speed  of  sound  water  and  air  respectively.  Once  in 
the  ocean,  this  wave  propagates  as  a  regular  plane  wave. 
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Under  certain  flight  conditions,  ray  paths  can  converge  to  form  superbooms.  Two  types 
of  focusing  can  lead  to  superbooms.  The  first  is  a  simple  focus  corresponding  to  a 
caustic.  In  three  dimensions,  a  caustic  is  a  two  dimensional  plane  which  intersects  the  air- 
water  interface  (or  ground)  in  a  line.  Typical  conditions  leading  to  superboom  caustic 
focusing  include  acceleration  and  steady  turns  (Figure  3-17).  Caustic  peak  pressure  focus 
factors  of  between  2  and  5  have  been  predicted  and  observed.  Focal  boom  footprints 
depend  on  the  maneuver  leading  to  the  caustic  but  are  generally  confined  around  the 
ground  track  below  where  the  maneuver  occurs.  The  second  level  of  focusing  known  as 
superfocusing  corresponds  to  a  cusp  between  two  smooth  caustics.  Superfoci  occur 
during  transient  maneuvers  such  as  a  turn  entry.  In  three  dimensions,  a  superfocus  cusp 
is  a  line  which  intersects  the  ground  at  a  point.  Superfoci  have  been  observed  with  shock 
focus  factors  approaching  10  times  observed  steady  level  flight  N-wave  pressures  at 
similar  speeds  and  altitudes.  The  extent  of  such  superfoci  is  limited  to  a  few  hundred  feet 
in  size. 


Flight  trad: 


Figure  3-17  Focusing  Due  to  Acceleration  (Maglieri  &  Plotkin,  1995) 
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Figure  3-18  Idealized  Inverted  Cycloid  Ocean  Wave 

All  analysis  to  this  point  has  assumed  a  flat  ocean  interface.  In  the  presence  of  a  rough 
interface,  wave  slope  may  allow  boom  energy  to  enter  the  ocean  at  angles  less  than 
critical.  Ocean  waves  can  be  characterized  as  inverted  cycloids  with  crest  to  crest  length 
L  and  crest  to  trough  height  H  relative  to  the  still  water  level.  The  wave  slope  can  be 
defined  as  angle  A  shown  in  Figure  3-18.  If  the  difference  of  the  Mach  angle  a  =  sin' 
*(1/M)  and  wave  slope  angle  A  is  less  than  the  critical  angle  0=13.3°  then  boom  energy 
may  be  injected  into  the  ocean  at  angles  less  than  critical. 

Bowditch  (1962)  provides  tables  which  relate  Beaufort  sea  state  number  and  fetch  in 
nautical  miles  to  wave  height  (H  in  feet)  and  period  (P  in  seconds)  for  a  fully  developed 
sea.  Bowditch  also  relates  the  wave  length  (L  in  feet)  to  the  period  P  as  :  L  =  5.12P  . 

Using  the  Bowditch  (1962)  wave  tables  and  the  relationship  between  P  and  L,  estimates 
for  angle  A  as  a  function  of  sea  state  and  fetch  can  be  computed.  Table  3-3  gives 
estimated  wave  slopes  for  fetch  distances  of  10  and  100  nautical  miles  for  Beaufort 
numbers  between  3  and  10.  For  reference,  Beaufort  number  3  corresponds  to  a  gentle 
breeze,  7-10  knots,  while  number  10  is  described  as  a  whole  gale  with  winds  in  the  48-55 
knot  range. 


Beaufort 

Number 

Wind  Speed 
(km/h) 

Wave  Angle  (A)  in 
degrees.  10  nmi.  Fetch 

Wave  Angle  (A)  in 
degrees.  100  nmi.  Fetch 

3 

12.9  - 19.3 

9.1 

2.8 

4 

20.9  -  29.0 

10.0 

4.6 

5 

30.6  -  38.6 

9.9 

5.8 

6 

40.2  -  50.0 

11.5 

7.1 

7 

51.5-61.2 

11.5 

8.2 

8 

62.8  -  74.0 

10.6 

9.5 

9 

75.6  -  86.9 

10.5 

10.2 

10 

88.5-101.4 

12.5 

10.8 

Table  3-3  Wave  Angle  Estimates  as  a  Function  of  Sea  State  and  Fetch 
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Boom  energy  will  propagate  into  the  ocean  if  the  wave  slope  angle.  A,  is  greater  than 
angle  B  =  a  - 13.3  where  a  is  the  Mach  angle  and  13.3  is  the  critical  angle.  Table  3-4 
shows  angle  B  as  a  function  of  Mach  number  and  shows  the  sea  states  where  the  wave 
slope,  A,  is  likely  to  be  larger  than  B  thus  allowing  boom  energy  to  propagate  into  the 
ocean  at  angles  below  critical. 


Mach  Number 

Angle  B  (degrees) 

=  a-13.3° 

Beaufort  Number 
10  nmi  Fetch 

Beaufort  Number 
100  nmi  Fetch 

3.5 

3.3 

>2 

>3 

3.0 

6.2 

>2 

>5 

2.5 

10.3 

>3 

>9 

2.0 

16.7 

>10 

>10 

1.5 

28.5 

>10 

>10 

Table  3-4  Speeds  and  Sea  States  Resulting  in  Boom  Waves  Below  Critical 


Table  3-4  indicates  that  for  speeds  below  Mach  2,  ocean  slopes  leading  to  propagating 
boom  injection  are  unlikely  for  all  but  the  most  severe  sea  states.  On  the  other  hand, 
wave  slopes  resulting  from  moderate  seas  are  likely  to  lead  to  sub-critical  boom  angles 
for  aircraft  speeds  greater  than  Mach  3.  This  analysis  is  supported  by  Sohn  et  al.,  2000 
who  conclude  “the  scattered  boom  signal  is  expected  to  be  negligible  until  vehicle  speeds 
reach  Mach  ~3.  At  Mach  ~3  very  rough  sea  states  have  the  potential  to  scatter  significant 
amounts  of  boom  energy  into  the  water  column.  Between  Mach  3  and  Mach  4.4  the 
magnitude  of  the  scattered  signal  will  increase  with  vehicle  speed  and  sea  state.” 

The  recent  measurements  of  Sohn  et  al.  (2000)  indicate  minimal  impact  of  rough  surface 
effects,  and  thus  support  the  theory,  but  only  for  the  case  of  low  sea  states  and  low  Mach 
numbers. 

As  noted  earlier,  waves  propagating  from  the  air  to  the  water  lose  about  30  dB  of 
intensity  due  to  the  impedance  mismatch.  Sub-critical  propagating  waves  may  also  suffer 
additional  losses  if  the  wavelengths  of  the  incident  energy  are  large  compared  to  the 
length  of  the  ocean  waves.  Sonic  boom  energy  peaks  are  in  the  2-20  Hz  range  which 
corresponds  to  wavelengths  in  the  17  - 172  meter  range.  Based  on  the  tables  from 
Bowditch,  1962,  typical  surface  wavelengths  are  on  the  order  of  10-30  m  for  low 
Beaufort  numbers  (3-4)  and  25-75  m  for  high  Beaufort  numbers  (8-9)  assuming  10-100 
nmi  fetches.  By  comparing  the  acoustic  wavelengths  associated  with  the  boom  energy  to 
those  expected  for  the  ocean  waves  it  is  apparent  that  rough  ocean  coupling  losses  may 
be  experienced  for  moderate  sea  states.  Sub-critical  wave  transmission  for  a  rough  ocean 
is  detailed  in  the  companion  subsonic  aircraft  noise  report  (Eller  and  Cavanagh,  2000). 

3.5  Methods  for  Estimating  Noise  Properties 

“The  United  States  Air  Force  is  required  to  analyze  the  environmental  impact  from  sonic 
booms  produced  by  supersonic  maneuvering  aircraft.  To  meet  this,  a  single  event 
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prediction  model  PCBOOM3  has  been  developed.  This  model  uses  full  ray  tracing 
propagation  theory  and  includes  the  effects  of  aircraft  maneuvers,  non-standard 
atmospheric  profiles,  and  winds.  Calculations  include  focused  superbooms  when  they 
occur.  Input  parameters  are  flexible  to  allow  general  supersonic  maneuvers.  The 
program  accepts  user  defined  F-functions,  and  has  F  functions  built  in  for  current  USAF 
aircraft.  The  model  predicts  the  location  and  waveforms  of  sonic  booms  on  the  ground. 
Predicted  sonic  boom  footprints  are  in  the  form  of  contours  (psf,  CSEL,  etc)  and 
isopemps,  or  as  tabulated  values  at  specific  points.”  (Plotkin,  1994).  The  aircraft 
trajectory  is  specified  by  an  initial  location,  altitude,  heading,  climb  angle  and  first  and 
second  derivatives  of  the  Mach  number,  heading  angle,  and  climb  angle.  Subsequent 
positions  are  projected  forward  in  time  using  constant  second  derivatives  with  an  option 
to  change  the  second  derivative  with  each  time  advance.  The  atmosphere  is  specified  by 
temperature  as  a  function  of  altitude  and  wind  speed  as  a  function  of  altitude  and 
direction.  The  temperature  profile  is  converted  into  a  sound  speed  profile  to  which  the 
wind  speed  profile  is  added.  This  method  of  including  wind  effects  is  a  standard 
technique,  but  treats  paths  perpendicular  to  the  ground  track  poorly. 

PCBoom3  is  used  for  single  sonic  boom  predictions.  “For  multiple  boom  environments, 
we  have  two  models,  CorBoom  and  BoomMap3.  CorBoom  predicts  the  yearly  average 
boom  levels,  Lcdn,  for  supersonic  operations  along  a  corridor.  BoomMap3  predicts  the 
Lcdn  contours  for  supersonic  operations  in  air  combat  maneuvering  instrument  areas 
(ACMI).  This  program  uses  the  ACMI  tracking  data  to  predict  the  sonic  boom 
environment.  Both  of  these  programs  used  a  simplified  ray  tracing  algorithm  developed 
by  Dr.  Kenneth  Plotkin  of  Wyle  Laboratories.”  (FICAN,  1994). 
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4.0  EXAMPLES  AND  CONCLUSIONS 


4.1  Examples 

Peak  pressure  and  N-wave  duration  can  be  estimated  using  the  formulas  in  equations  [1] 
or  computed  using  PCBoom3.  Equations  [1]  require  aircraft  length  and  mass  which  are 
summarized  in  table  4-1  below  for  several  typical  aircraft. 


Type 

Mass  (klbs) 

Length  (ft) 

B-l 

Variable  Sweep 

450 

147 

F-4 

Large  Fighter 

45 

60 

F-14 

Variable  Sweep 

55 

62 

F-15 

Fixed  Wing  Fighter 

65 

64 

F-16 

Fixed  Wing  Fighter 

35 

49 

F-18 

Fixed  Wing  Fighter 

39.4 

56 

F-22 

Fixed  Wing  Fighter 

48 

67 

F-104 

Small  Fighter 

26.5 

55 

F-lll 

Variable  Sweep 

76.5 

74 

T-38 

Small  Fighter 

20 

47 

Concorde 

Supersonic  Transport 

387 

190 

Shuttle 

Blunt  Lifting  Body 

187 

121 

Table  4-1  Selected  Aircraft  Mass  and  Lengths  (PCBoom3  &  WWW) 


Equations  [1]  also  require  atmospheric  pressure,  Pv,  at  the  aircraft  altitude  h.  Pv  may  be 
crudely  estimated  using  the  exponential  form  presented  in  [1]  or  may  be  taken  from  more 
precise  measurements.  The  examples  below  use  pressures  from  the  standard  U.S. 
atmosphere  as  summarized  in  Table  4-2  (U.S.  Standard  Atmosphere,  1976). 


Altitude  (km) 

Pressure  (atm) 

Pressure  (Pa) 

0 

1.0 

101,327 

5 

0.533 

54,021 

10 

0.261 

26,437 

15 

0.119 

12,045 

20 

0.054 

5,475 

25 

0.025 

2,511 

Table  4-2  Standard  Atmospheric  Pressure 


Two  cases  are  presented  below.  The  first  involves  a  450  klb.,  147  foot  long,  B-l  variable 
sweep  wing  bomber  flying  at  Mach  1.2.  Assume  the  aircraft  is  in  steady  level  flight  over 
a  flat  ocean.  The  peak  pressures  and  N-wave  durations  as  estimated  by  Carlson’s 
equations  [1]  and  PCBoom3  are  shown  in  Figure  4-1  as  a  function  of  aircraft  altitude. 

The  PCBoom3  pressure  results  were  scaled  by  a  factor  of  2.0/1. 9  to  account  for  the 
surface  reflection  factor  expected  over  water  (2.0)  as  opposed  to  average  ground  levels 
(1.9).  In  general  the  two  models  agree,  with  Carlson  predicting  slightly  higher  peak 
pressures,  especially  at  low  altitudes.  Figure  4-2  shows  peak  pressure  and  N-wave 
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duration  variability  as  a  function  of  Mach  number  for  a  B-l  bomber  in  steady  level  flight 
at  10  km  altitude.  The  peak  pressure  increases  about  25%  as  the  aircraft  speed  increases 
from  M=1.2  to  M=3.0.  Equations  [2]  can  be  used  to  estimate  peak  pressures  as  a  function 
of  depth  in  the  ocean.  As  mentioned  earlier,  the  lateral  spread  of  the  primary  boom 
carpet  can  be  estimated  by  equating  the  lateral  spread  in  miles  to  the  aircraft  altitude  in 
thousands  of  feet. 

For  a  specific  example  of  a  flyover  time  signature,  consider  a  B-l  bomber  traveling  at 
Mach  1.2  at  an  altitude  of  10  km.  Based  on  figure  4-1,  the  N-wave  peak  pressure  at  the 
air-water  interface  is  3.7  lb/ft2  =164.8  dB  which  includes  the  surface  reflection  factor. 
The  corresponding  N-wave  duration  is  0.325  seconds  which  for  an  aircraft  speed  of 
1.2*333m/s  leads  to  an  N-wave  length  of  129.9  m.  Using  these  values,  the  sound 
pressure  as  a  function  of  time  from  the  center  of  the  N-wave  and  depth  is  shown  in  figure 
4-3. 


The  second  case  study  involves  a  26.5  klb  ,  55  foot  long,  F-104  small  fixed  wing  fighter 
flying  at  Mach  2.0.  Figure  4-4  shows  peak  pressure  and  N-wave  duration  for  a  F-104  in 
steady  level  flight  at  Mach  2.0  as  a  function  of  altitude.  The  Carlson  and  PCBoom 
models  agree  very  well,  with  the  Carlson  peak  pressures  being  slightly  higher  at  low 
altitudes.  Figure  4-5  shows  peak  pressure  and  N-wave  duration  for  a  F-104  in  steady 
level  flight  at  10  km  altitude  as  a  function  of  Mach  number.  As  with  the  B-l  bomber  case 
above,  the  Carlson  N-wave  durations  fall  slightly  below  the  PCBoom3  estimates.  Figure 
4-6  shows  a  10  km  altitude,  Mach  2.0  F-104  flyover  time  signature. 

Table  4-1  provides  N-wave  duration,  peak  pressure  and  energy  flux  density  level  at  the 
ocean  surface,  50  m,  and  100  m  below  the  surface  for  a  variety  of  aircraft  flying  at 
several  altitudes  and  speeds.  The  energy  density  levels,  peak  pressures,  and  N-wave 
durations  were  computed  using  the  work  of  Carlson,  1978  (equation  [1]  and  figures  8-9) 
and  Cook,  1970  (equations  [2]).  Energy  density  levels  were  computed  by  numerical 
integration  where  the  pressure  time  series  was  sampled  at  Nyquist  frequency  of  at  least 
1500  Hz.  As  noted  in  section  3.2,  spectrum  levels  are  at  least  40  dB  down  from  their 
peak  above  1000  Hz  for  surface  N-waves.  For  sonic  boom  arrivals  sampled  below  the 
surface,  spectrum  levels  decay  faster  than  the  6  dB/octave  seen  in  the  surface  N-wave. 
Therefore,  spectrum  levels  below  the  surface  will  be  more  than  40  dB  down  from  their 
peak  values  above  1000  Hz.  With  this  information,  the  peak  and  energy  density  values 
presented  in  table  7  can  be  used  for  frequencies  less  than  1000  Hz.  Frequencies  greater 
than  1000  Hz  will  be  at  least  40  dB  lower  than  the  table  values. 
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Figure  4-1  B-l  Bomber  at  Mach  1.2 


Figure  4-2  B-l  Bomber  at  10  km 
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Table  4-1  Sonic  Boom  Levels  for  Selected  Aircraft 


4.2  Remarks  and  Conclusions 

At  the  most  fundamental  level,  consider  the  acoustic  principle  for  fluids  that  across  any 
plane  the  pressure  must  be  continuous  (as  must  be  the  normal  component  of  the  particle 
velocity).  This  principle  is  essential  for  estimating  acoustic  propagation  from  air  to 
water  and  vice  versa,  and  is  invoked  in  virtually  all  referenced  treatments  of  the  problem. 
From  a  practical  point  of  view,  it  means  that  twice  the  pressure  in  air  at  the  air-sea 
interface  bounds  the  pressure  in  the  water  just  below  it.  Hence  a  sonic  boom  of  peak 
pressure  10  psf  at  the  ocean  surface  becomes  an  impulsive  wave  in  water  with  a 
maximum  peak  pressure  of  20  psf. 

The  issue  then  is  whether  the  impulsive  wave  in  the  water  propagates  in  depth  and  range, 
or  not.  In  most  cases  of  practical  supersonic  flight,  the  wave  is  not  a  propagating  wave, 
but  rather  a  wave  that  decays  in  pressure  relatively  rapidly  in  depth  or  range.  As 
reviewed  in  Section  3,  much  research  has  gone  into  this  problem,  including  recent  work 
by  Sparrow  (1995),  Rochat  and  Sparrow  (1995),  Cheng  et  al.  (1996,  1997),  and  Sohn  et 
al.  (2000),  and  others. 

In  some  very  special  cases  (e.g.,  Mach  number  greater  than  4.5,  special  aircraft 
maneuvers  that  momentarily  steepen  the  shock  wave  angle  at  the  ocean  surface,  perhaps 
the  case  of  a  very  rough  sea)  sonic  boom  energy  may  in  fact  become  propagating  energy 
in  water.  But  again,  for  practical  geometries,  twice  the  pressure  at  the  interface 
approximates  the  maximum  pressure  that  will  be  observed  in  range  or  depth. 

Now,  a  very  loud  sonic  boom  wave  at  the  ocean  surface  would  be  one  of  order  10  psf. 
Anything  above  50  psf  would  be  difficult  to  generate  (low-flying  missiles  having  been 
estimated  to  yield  peak  pressures  on  the  ground  as  great  as  35  psf).  Hence,  100  psf  is 
proposed  here  as  an  upper  bound  that  anticipates  the  worst  cases.  For  comparisons 
below,  note  that  100  psf  equals  about  0.7  psi,  or  about  194  dB  (re  1  |iPa). 

As  discussed  at  length  in  the  companion  volume  on  thresholds  for  injury  and  harassment 
of  marine  life  (Cavanagh  and  Laney,  2000),  the  most  commonly  used  thresholds  for 
harassment  of  marine  mammals  and  sea  turtles  by  impulsive  noise  are: 

12  psi  peak  pressure,  and/or 

182  dB  (re  1  (xPa2-s)  energy  flux  density  level  in  the  most  energetic  1/3  octave 
band  above  10  or  100  Hz  (depending  on  the  species). 

Thus,  even  a  worst  case  situation  for  sonic  booms  will  yield  a  peak  pressure  in  water  well 
below  the  12  psi  threshold.  As  for  the  energy  metric,  the  calculations  in  this  volume  of 
the  energy  of  sonic  boom  waves  in  water  yield  levels  of  order  180  dB  or  less  for  the  100 
psf  boom.  Table  4-1  can  help  to  put  this  is  perspective. 
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It  is  also  worthwhile  to  consider  the  inverse  problem.  To  approach  typical  harassment 
levels  for  marine  mammals,  we  must  have  about  12  psi  peak  pressure  in  the  water.  This 
means  that  peak  pressure  of  the  same  order  must  occur  at  the  surface,  in  air.  But  half  the 
pressure  in  air  is  6  psi,  or  about  900  psf.  Note  that  a  sonic  boom  peak  pressure  of  50  psf 
at  the  ocean  surface  is  considered  extraordinary. 

For  comparison  to  other  impulsive  noises  in  water,  consider  the  standard  Navy 
underwater  SUS  charge  (about  1.8  pounds  of  TNT  equivalent).  At  1  meter  from  the 
charge,  the  peak  pressure  is  about  260  dB  (re  1  microPa)  -  about  1000  psi  and  144,000 
psf.  As  far  away  as  1  km,  the  pressure  will  be  about  1  psi,  or  144  psf.  To  reach  the  level 
of  a  loud  sonic  boom  at  the  surface  (say  10  psf),  the  observer  would  have  to  be  about  7 
km  away.  Yet  SUS  charges  are  routinely  used  throughout  the  world’s  ocean  and 
assessed  by  Navy,  with  agreement  of  the  regulators,  as  an  insignificant  risk  to  protected 
marine  species. 

It  is  evident,  then,  that  by  today’s  standards  sonic  booms  present  essentially  no  risk  of 
harassment  for  protected  marine  species  in  water. 

Nonetheless,  estimates  of  the  type  made  in  this  volume  are  essential  justifications  for 
arguments  of  insignificant  impact  (as  would  be  used  in  compliance  documents).  They 
also  anticipate  the  event  that  (as  happens  every  few  years)  the  favored  harassment  and 
injury  thresholds  are  changed  (from  the  view  of  the  regulators).  As  recently  as  four 
years  ago,  mammal  harassment  thresholds  used  in  Navy  ship  shock  trials  were  stated  to 
be  as  low  as  160  dB  (energy  level).  As  Table  4-1  indicates,  such  levels  can  well  occur 
for  selected  cases.  Detailed  risk  analyses  for  specific  animal  populations  and  locations 
would  then  have  to  be  developed,  and  mitigation  schemes  designed  to  reduce  risk. 
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